A statistical model for the 7-ray variability of the Crab nebula 
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ABSTRACT 

A statistical scenario is proposed to explain the 7-ray variability and flares 
of the Crab nebula, which were observed recently by the Fermi/LAT. In this 
scenario electrons are accelerated in a series of knots, whose sizes follow a power- 
law distribution. These knots presumably move outwards from the pulsar and 
have a distribution in the Doppler boost factor. The maximal electron energy is 
assumed to be proportional to the size of the knot. Fluctuations at the highest 
energy end of the overall electron distribution will result in variable 7-ray emission 
via the synchrotron process in the ~ 100 MeV range. Since highly boosted larger 
knots are rarer than smaller knots, the model predicts that the variability of 
the synchrotron emission increases with the photon energy. We realize such a 
scenario with a Monte-Carlo simulation and find that the model can reproduce 
both the two 7-ray flares over a period of ~ year and the monthly scale 7-ray 
flux fluctuations as observed by the Fermi/LAT. The observed 7-ray spectra in 
both the steady and flaring states are also well reproduced. 

Subject headings: radiation mechanism: non-thermal — pulsar wind nebula: 
individual: Crab — gamma rays: theory 



1. Introduction 



The Crab pulsar wind nebula is powered by its central pulsar born from a supernova 
explosion in 1054. It is a very luminous source in almost all wavelengths, from radio to the 
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very high energy (VHE) 7-ray bands. The broadband non-th ermal emission spectrum is 



well modeled by a synchrotron-inverse Compton (IC) scenario flAtoyan fc Aharonianlll996 



Meyer et al.ll2010l ) with the transition from the synchrotron to IC component occurring at a 
few hundred MeV flAbdo et al.ll2010bl ). The overall emission seems to be steady, so that the 
Crab nebula has been adopted as a standard candle in high energy astrophysics to calibrate 
observations from different instruments. 

However, detailed images of the Crab nebula in optical and X-ray bands indicated 
dynamical structures at small scales. Observations by the Hubble Space Telescope revealed 
wisps and knots in the nebula, and resolved some substructures of the highly variable wisps 
f lHester et al.l Il995[ ). X-ray o bservations by ROSAT and Chandra u ncovered a jet-torus 

and the equatorial 



I2OOO 



structure of the inner nebula ( iHester et al.lll995t IWeisskopf et al. 
ring is found moving outwards with a speed of ~ 0.5c f lHester et al.ll2002[ ). At higher energies 



in the 7-ray band, detailed structures of the nebula generally can not be resolved. However, 
COMPTEL and EGRET observations in dicated that the synchrotron radiation in 1 — 150 
MeV is variable on a timescale of ~ 1 yr ( iMuch et al.lll995l : Ide Jager et al.lll996l ). 



In September 2010, AGILE collaboration reported a 7-ray flare above 100 MeV from 
the direction of the Crab nebula, which lasted for about 3 days (ATel #2855; Tavani et al. 
2011). The flux during the flare period is about 2 — 3 times higher than the average one. 
This flare was soon confirmed by the Fermi/LAT collaboration (ATel #2861), and archiv e 
search of the Fermi/LAT data revealed another flare in February 2009 (lAbdo et al.ll2010al ). 
There were many other simultaneous or follow-up measurements for the flare in September 
2010 in X-ray, optical, infrared and radio bands. However, no significant flux enhancement 
was discovered. In the VHE 7-ray energies, ARGO-YBJ collaboration claimed the detection 
of a flux enhancement around TeV, with a possibly longer duration (ATel #2921). However, 
MAGIC and VERITAS observed the source with a shorter duration during the flare phase 
and did not reveal any enhancement in flux (ATel #2967, #2968). 

It was proposed that the flares was due to synch rotron emission of ultra- relativistic 
electrons with energies up to ~PeV flAbdo et al.ll2010al ). Considering the fact that the rest 
frame synchrotron radiation in a magnetic field dominated acceleration regime can not exceed 
_70 MeV due to the fast synchrotron cooling of high energy electrons ( iBlumenthal fc Gould 



I970I : iLyutikovl I2OIOI ) , the observed GeV emission implies Doppler s hift of the radiatio n 
region or other acceleration rn echanisms instead of shock acceleration ( lAbdo et al.ll2010al ). 
Komissarov fc Lyutikovl (120101 ) proposed that the 7-ray va riability (flare ) origi nated from the 
"inner knot" of the Crab nebula ("knot 1" as defined in iHester et al.l (119951 )). with mildly 
Doppler-boosted emission. The instability of the termination sh ock may cause the vari- 
ability as revealed by magnetohydrodynamic (MHD) simulations. iBednarek fc Idea (120 lOf ) 
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suggested that electrons are accelerated in a region behind the shock, and the variability 
was attributed to changes in the maximum energy of accelerated electrons, electron spectral 
index or the magnetic field. 

In this Letter we employ a statistical approach to model the 7-ray variability of the 
nebula. Fermi/LAT observations have shown that the low-energy synchrotron component 
i s variable on mont hly time scale, while the high-energy IC component seems to be stable 



( lAbdo et al.ll2010al ). Since the synchrotron 7-rays are produced by the highest energy elec- 
trons, these observations indicate that fiuctuations at the high-energy end of the electron 
distribution might be responsible for the variability. It is natural to expect that events that 
can generate the highest energy electrons are rarer, and therefore would lead to the largest 
fiuctuation. Therefore the variability and fiares in the sub-GeV 7-rays can be simply due 
to the statistical fiuctuation of the highest energy electrons achievable in the electron accel- 
erators. Lower energy electrons do not suffer from significant fiuctuations since many more 
accelerators can contribute to them simultaneously, which gives rise to a "steady-state" emis- 
sion in both the lower energy synchrotron component and the higher energy IC component. 



2. The Model and Results 

We build a physically possible model to investigate the above picture in more details. We 
assume that electrons are accelerated in a series of knots with a power-law energy spectrum 
Fi[E) oc with the maximum electron energy -E^ax being proportional to the size of the 

ith knoiQ. We further assume that the size of knots has a power-law distribution P(rj) oc . 
Adding the contribution from all these knots together, noticing that the normalization of 
the synchrotron spectrum of the knots scales as rf, we get the total electron spectrum as 

F{E) oc J F.(E)e(ELx -E)x rfP{r,)dn 

(X J P-"^e(i?Lx -E)x ElJ-^dEl,^ 

= E I El^^^ dE^a^, (1) 
Je 



^The maximum energy of electrons will reach a critical value when the synchrotron cooling is essential. The 
critical Larmor radius is estimated as i?£ ~ \l "eB — 2 x 10~'^(i3/mG)~"'^'^pc, which just corresponds 
to the size of the emission region as indicated by the day-scale flare. Therefore we assume the largest size of 
the knots is i?^, and for all other knots smaller than i?^ the scaling relation between -E^^x ^-^id can well 
hold. 
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where Q{x) is the Heaviside step function, -Emax is the maximum energy of the largest knot. 
If /3 > 4, the above integral of the last step is approximately proportinal to given the 
maximum energy of the largest knot -Emax ^ -E; if /3 < 4 the integral is nearly a constant 
determined by -Emax- Therefore the electron spectrum is a power law F{E) oc with 
power-law index a = max(Q!e + /3 — 4, ae) and a cutoff determined by the largest knot. 

The maximal energy of synchrotron emission in the magnetic field dominated accelera- 



tion regime is ~ 70 MeV (IBlumenthal &: Gouldlll970l ) (defined by equating the synchrotron 



cooling time scale to the gyro-period, the minimum of the acceleration time scale). For 
smaller accelerators, the maximum energy is also limited by the size of the accelerator. In 
order to account for the observations in the ~ 100 MeV range, we therefore further employ 
a mild Doppler boost factor. We assume a Gaussian distributed Lorentz factor of all the 
knots with a mean value F = 2.0 and a standard deviation cr = 0.25. Such a Lorentz fac- 



tor is consistent with the upper limit of the typical velocity of the jet (lAbdo et al.ll2010af ). 
We assume that the angle 6 between the knot motion and the line-of-sight is randomly dis- 
tributed. The Doppler factor 6 = 1/F(1 —v cos6) (where v is the velocity of the knot in unit 
of c) therefore is distributed within the range 0.18 — 5.5 for a 3a distribution of F. Only 
knots with large enough size and large enough Doppler factor can contribute to high energy 
synchrotron radiation to account for the flares observed by the Fermi/LAT. 

We realize such a picture by a Monte-Carlo simulation. We directly simulate the syn- 
chrotron spectrum in each knot instead of the electron spectrum. The synchrotron spec- 
trum of each knot in its co-moving frame is adopted as z/'Fj,/ cx z/'""" exp(— z/'/z/^g^), where 
'^max (-^max)^ ^ OLy = ipi — 3)/2 is the syuchrotrou spectral index, where a is the 
power-law index of the overall electron distribution. The magnetic field in each knot is not 
explicitly used in the calculation. The magnetic field will affect the cooling break of the low 
energy synchrotron spectrum as well as the high energy cooling time scale. To keep the basic 
framework of the present study, we may expect the magnetic field not to vary significantly 
with respect to the knot size. 

The energy loss of electrons in the knots can be important and needs to be consid- 
ered. For instantaneous injection of accelerated electrons, the probability of observing 
electrons at a given energy is proportional to the corresponding synchrotron cooling time 
tc oc . For low energy electrons with cooling time longer than the age of the Crab nebula 
tage ~ 10^ yr, these electrons will survive and accumulate. Therefore we can construct a 
probability function P(-E') = 1/(1 + tageAc), describing the effect of cooling on the electron 
spectrum. The cooling time of electrons producing synchrotron photon with energy e is 
tc — 1.5(S/mG)~^'^(e/keV)~°'^5"°'^ yr. For an average magnetic field B ^ 0.1 niG, a mildly 
Doppler factor 5 ~ 1, electrons producing photons with an energy e ~ 2.5 eV have a cooling 
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time comparable to tagc- The observed spectrum of Crab nebula has a break near 2.5 eV 
(see below Fig. [1]). 

In this work we mainly focus on the high energy part of the spectra, e.g., from X-ray 
to 7-ray band, then P{E) oc E~^. Assuming that the accelerated electrons have a spectrum 
index the electron spectrum index in the synchrotron cooling dominated regime is then 
ae = a"-' + 1. Here we adopt ap = 2.6 wh ich corresponds to a™-" = 1.6, as required by fitting 



the radio-optical data (IMeyer et al.ll2010l ). The maximum value of t^'^^^ is taken as 55 MeV 



in the co-moving frame (corresponding to the largest knot). This value is derive d according 



to a fi t to the average Fermi/LAT spectrum of the synchrotron component (lAbdo et al. 



2010bl ). given the Lorentz factor distribution. The power-law index of the size distribution 
is adopted as /3 ~ 4.8. Adding all the contribution from the knots together we can get a 
total synchrotron spectrum z/'Fj,/ oc z/'~°'^, which can reproduce the observed optical- MeV 
band data of the Crab nebula (Fig. [1]). Note that these parameters, i.e., F, i^^ax' '^e and (3 
are not uniquely determined. What we adopt in this work is an illustration of the model. 
Once the data from optical to 7-ray band are reproduced, the discussion below is basically 
unchanged if a different set of parameters is taken. 

For each knot we randomly generate an angle 6 and a Lorentz factor F based on the 
assumed Gaussian distribution, and calculate the Doppler factor 6. Then the frequency of 
synchrotron emission from each knot is shifted from u' to Su', and the fiux is enhanced to 
6^F^r. The contribution from all the knots are summed up to obtain the total emission 
spectrum. 

A series of simulated synchrotron spectra are shown in Fig. [H which are compared with 
the observational data of the Crab nebula. Note that for the INTEGRAL data we adopt a 
fiux nor malization fac t or 0.6 in order to be consistent with the extrapolation of COMPTEL 



data. In iMeyer et al.l (j20ld ) a 0.78 factor was adopted to match the XMM- Newton data. 
The difference of flux normalization may be due to discrepancies in calibration of different 
detectors. The fact that the INTEGRAL data may contain emission from the Crab pulsar 
may also contribute to the difference. For comparison we show the fltting results to the broad 
band data with the blue thick line, with a broken power-law electron sp ectrum. This result is 



(1996 


) and 


Mever et al. 


(2010) 



small fluctuations in the X-ray energy band, but fluctuate signiflcantly at higher energies 
(> 1 MeV). The average of the simulated 7-ray flux would be consistent with the long term 
Fermi/LAT and COMPTEL data. The fluctuations can be responsible for the observed 
variabilities and flares. 



The very large flare at September 2010 that is not well reproduced by these simulations 
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Fig. 1. — Spectral energy distribution of the Crab nebula. The thin green lines de- 
note the simulated synchrotron spectra of 30 realizations. The thick blue line denote 
the fit to the multi-wavelength steady state observational data, which can be under- 
stood as the average result of many realizations in our model. The black-dashed curve 
denotes the contribution from a maximum size knot with slightly specific parameters 
(see the text), which is shown to be able to explain the larg e flare in Septembe r 2010 . 
The reference s of the observational da ta are: radio-o ptical f|Macias-Perez et al.l |2010| ). 
INTEGRAL J.Tourdain fc Ronuesl boOflh. COMPTEL JKuiper et all boOlh Fem ii/LAT 
flAbdo et allboiObl ). MAGIC f lAlbert et al.lboosh . HEGR A JAharonian etil] 120041 ). HESS 



(Aharonian et al. 



2006 ). Tibet-AS7 ( Amenomori et al. 20091 ). and Fermi/LAT flares 



(lAbdo et al.ll2010af ). 



can be regarded as an event with very small probability. Such an event may be due to a large 
knot with specific parameters different from what used in the model. The black-dashed line 
in Fig. [T] shows a possible reproduction of the large flare, produced by a knot with maximum 
size, Doppler factor 6 = 5.5, cutoff energy of synchrotron radiation in its comoving system 
70 MeV and an additional normalization factor 0.2 of the flux. Detailed modeling of the 
knot formation and its evolution is expected to better address this issue. 

We pick out two realizations and show the synchrotron sky-maps at 10 keV (left) and 
100 MeV (right) in Fig. |2l respectively. To generate a sky-map, we designate a.n x — y plane 
with 50 X 50 pixels. We randomly assign a position for each knot in the x — y plane, with 
the largest knot having a size of 5 x 5 pixels. We then sum up contributions of all the knots 
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located in each pixel and obtain the sky-map. It is shown that in the X-ray band (left panels) 
the skymaps are relatively smooth with some hot spots due to several large knots. These hot 
spots may correspond to the Chandra and Hubble knots of the Crab nebula (ATel #2882, 
#2903). The total flux of X-ray emission, however, is statistically steady from time to time. 
The width of the total X-ray flux distribution at ~ 10 keV is estimated to be of the order 
1%, which is consistent with the upper limit of 5.5% obtained by Swift/BAT (ATel #2858). 
The ~ 100 MeV 7-ray sky-map could differ significantly from each other. In some cases, 
one or several large knots would dominate the total flux. The predicted energy-dependent 
spatial scales due to fluctuations may not be ready to test with the available high energy 
telescopes, but higher angular resolution observations in the future will help. 

According to our model, the IC component should not show significant variability except 
in the high energy regime (> 100 TeV), depending on the maximal electron energies. For 
example, for a magnetic field B 1 mG, the maximal Lorentz factor of electron can be 
as high as 7max ~ ^/Qire/arB ^ 4 x 10^. Therefore only the IC photons with energies 
higher than several hun dred TeV would have re latively large fluctuations. This conclusion 



is different from that of iBednarek fc Idea (120101 ). who suggested that the IC 7-ray spectra 



also vary significantly at multi-TeV energies. 

To investigate fluctuations of emission in detail, we show in Fig. [3]the light curves of the 
simulated synchrotron fluxes in the ~ 100 MeV energy band. Each bin of the light curves 
represents an independent realization, and the length of the time bin is proportional to the 
number of knots. We normalize the time scale so that over a year there is an enhancement in 
flux by at least a factor of 5-6, which may be responsible for a flare, and plot the weekly and 
monthly bin lightcurves in the upper and lower panels of Fig. [3l respectiyely. Fo r comparison 



we also plot the Fermi/LAT observed monthly light curve ( lAbdo et al.ll2010al ) in the lower 
panel of Fig. |3l It can be seen that the predicted fluctuation level of the simulated results is 
very similar to that observed by Fermi/LAT. In order to show this more clearly, we plot the 
histogram of the monthly bin flux distribution for the simulated results along with that of 
Fermi/LAT observation in the left panel of Fig. HI These two distributions are indeed similar. 
In the right panel of Fig. H] the histogram of flux distribution for weekly bin is also plotted. 
A concentration toward low-flux is evident, which is distinct from a symmetric distribution 
around the average value. This is because in the short time scale limit, the distribution of 
the gamma-ray fluxes is determined by the distributions of the knots (power-law) and the 
Doppler factors. Only the knots with both large enough sizes and high enough Doppler 
factors can produce large gamma-ray fluxes. Therefore the probability to give large fluxes 
is very small, and the flux distribution concentrates towards low-flux values. In the limit 
of long time scale, the distribution should be Gaussian as a result of central limit theorem. 
Such a prediction can be tested by a detailed analysis of the Fermi/LAT data. 
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Fig. 2. — Upper and bottom panels represent two realizations of the synchrotron sky-maps at 
10 keV (left) and 100 MeV (right) energies. The upper one shows a realization of a possible 
7-ray flare, while the lower one shows a realization with relatively weak 7-ray emission. 

3. Conclusion and discussion 

In summary, we propose a statistical picture to explain the observed 7-ray variability and 
flares of the Crab nebula, using the fluctuations of the highest end of the electron spectra. 
The electrons are thought to be accelerated in a series of knots, with a size distribution 
P(rj) cx r~'^ and a distribution of the Doppler factor. The maximal energy of the electrons 
in the co-moving frame is assumed to be proportional to the size of the knots. Thus the 
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Fig. 3. — Light curves of the simulated synchrotron flux above 100 MeV. The upper panel 
shows the results with weekly bin and the bottom panel for monthly bin. The vertical 
axes in each panel indicate the average v alue of the simulat ions, which is normalized to the 
observational value ~ 6.2x10"^ cm~^ s~^ f Abdo et ahlboiOa ). The Fermi/LAT observational 
data with monthly bin are also plotted in the lower panel. 



rare knots with large sizes and high Doppler boosts may generate electrons up to ~PeV, and 
hence, be responsible for the observed 7-rays flares at Fermi/LAT. On the other hand, the 
low energy electrons are generated by many smaller knots. The average effect smooths the 
fluctuation, so that the low energy synchrotron component and the IC component do not 
change significantly. This scenario can naturally explain both the variability in the MeV-GeV 
band, and the relatively steady emission in lower (optical, X-ray) and higher (TeV) energies. 
The expected variability of the monthly bin fluxes above 100 MeV are well consistent with 
that observed by Fermi/LAT. The two large 7-ray flares can also be naturally accounted for 
without additional assumptions. 

In the simulation, we do not consider the detailed structure (e.g. jet, torus) of the 
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Fig. 4. — Histograms of the synchrotron flux distribution above 100 MeV, for monthly bin 
(left) and weekly bin (right), respectively. In the left panel the Fermi/LAT observation result 
is also shown (dashed histogram). 

nebula. The knots are assumed to have a uniform and isotropic distribution. Considering a 
more realistic morphology of the nebula would not affect our conclusion noticeably, as long as 
the distributions of the knot sizes and Doppler factors are similar to those introduced here. In 
more general terms, these knots may be interpreted as individual electron acceleration events 
distributed throughout the nebula. We point out that at low energies contributions from 
large amount of knots are equivalent with previous studies that electrons are accelerated 
continuously inside the whole nebula. Although the model does not discuss the absolute 
time scales of flares and 7-ray fluctuations from the first principle, they may be related to 
timescales of developing MHD instabilities at different spatial scales. 
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